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Ames Research Center 


SUMMARY 


The complete flightpath optimization problem for minimum fuel consumption from 
takeoff to landing including the initial and final turns from and to the runway head- 
ing is solved. However, only the initial and final segments which contain the turns 
are treated, since the straight-line climbout, cruise, and descent problems have 
already been solved. The paths are derived by generating fields of extremals, using 
the necessary conditions of optimal control together with singular arcs and state 
constraints. Results show that the speed profiles for straight flight and turning 
flight are essentially identical except for the final horizontal accelerating or 
decelerating turns. The optimal turns require no abrupt maneuvers, and an approxima- 
tion of the optimal turns could be easily integrated with present straight-line climb- 
cruise-descent fuel-optimization algorithms. Climbout at the optimal IAS rather than 
the 250-knot terminal-area speed limit would save 36 lb of fuel for the 727-100 
aircraft . 


INTRODUCTION 


The escalating costs of fuel have increasingly focused attention on the problem 
of minimizing fuel consumption of commercial aircraft. Recent work in aircraft 
guidance has demonstrated that onboard optimization of aircraft trajectories offers 
an efficient method of fuel conservation. 

Aircraft trajectory-minimum fuel problems can be divided into two classes: 
en route problems, with trajectory lengths of 50 n. mi. and longer; and terminal-area 
problems with lengths of 3 to 50 n. mi. For the en route problem, an onboard algo- 
rithm for optimum climb-cruise-descent has been developed by applying optimal control 
theory (ref. 1). The solution was developed for the vertical plane only, since 
horizontal maneuvers are not a significant feature of en route flight. 

The general class of terminal-area trajectory problems are more difficult to 
solve, because vertical and horizontal maneuvers involving speed, altitude, and head- 
ing changes occur simultaneously and because they are of comparable significance in 
influencing fuel consumption. 

In this report we study the subclass of terminal-area climbout and descent tra- 
jectories which are part of an overall optimal trajectory that includes an en route 
straight-line optimal cruise segment. These maneuvers complete the optimal cruise 
trajectory problem of reference 1 by including the initial and final turns from and 
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to runway headings. For this subclass of terminal-area problems, essentially all 
possible extremals can be found, and the problem of global optimality encountered in 
the general class of terminal-area turning flightpaths can be resolved (ref. 2). 

The results presented in this paper show that for fuel-efficient turning climbout 
and descent guidance, the altitude-speed profile is nearly the same as that generated 
for straight-line flightpaths. In the special case of straight-line flightpaths, the 
results presented in this paper are identical to those presented in reference 1. 

Further research is needed on a second class of problems, that of optimal capture 
trajectories. These are trajectories which end by aligning the aircraft with the 
Instrument Landing System (ILS) localizer at a specific point and heading and begin at 
any other point and heading. These problems arise especially in the descent phase 
when the aircraft does not arrive in the terminal area at the proper state to complete 
an optimal approach or, for air-traffic-control reasons, is not permitted to complete 
the optimal maneuver. The approach in previous studies has been to synthesize and 
sometimes optimize vertical (ref. 1) and horizontal trajectories (refs. 3-6) sepa- 
rately and then combine them heuristically (refs. 7-9). In further research one can 
use extremals to study samples of optimal three-dimensional trajectories, which 
include speed changes, and thus verify and refine existing heuristic solutions. 
However, this problem is. not addressed in this paper. 


PROBLEM STATEMENT 


The problem is to minimize the fuel consumption over the complete flightpath 
from takeoff to touchdown, F*p. In order to solve this problem, the equations of 
motion for the vehicle are needed and the performance index must be formulated. 

The point-mass equations of motion for an aircraft in steady-state flight 
(weight = lift /cos <J>) , assuming small flightpath angle, constant mass, and coordi- 
nated turns are: 


= v cos 


(i) 

= v sin 


(2) 

= vy 


(3) 

= -gu/v 


(4) 

= g(T - 

D - Wy)/W 

(5) 


Here, x and y are the coordinates in the horizontal plane, h is the vertical coor- 
dinate, \p is the heading angle measured counterclockwise from the x-axis, v is the 
ground speed, g is the gravitational constant, W is the weight, and D is the drag. 
The control variables are: 


1. Thrust, T, where thrust is constrained to be between idle and maximum 
allowable thrust: 


T 


idle - 


T < T 
- max 
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2. Flightpath angle, y, which is constrained in accordance with present commer- 
cial airline practice: 


Y . 
mm 

VI 

>• 

VI 

Y 

max 


Y • 
mm 

= o°. 

Y =5° 

’max 

for climbout 

Y . 
mm 

= -4° 

, y =0° 

'max 

for descent 


3. The tangent of the bank angle, u, which is constrained in accordance with 
present practice for commercial autopilots. 

-u < u < u u = tan d> , <b = 30° 
m ~ ~ m m 

The only state constraint that was applied was on airspeed. FAA regulations state 
that the indicated airspeed in the terminal area must not exceed vj = 250 knots 
(422 ft /sec); when converted into true airspeed (which is equal to the ground speed 
in the assumed no-wind condition) this translates into a state variable constraint 
for the speed v: 

S = v - Vj(l + qh) <0; = ^22 ft/sec; q = 0.162xl0 -1 * (6) 

The next problem is to develop the performance index for the low-altitude turn- 
ing segment. This requires some preliminary discussion of various approximations. 

The intended solution for the complete minimal-fuel turning flightpath is to connect 
the optimal solution for the straight-line flightpath from reference 1 (cut off below 
an altitude of 10,000 ft) with an optimal turning climbout and descent (connected to 
the straight-line path at the 10,000-ft altitude). From Bellman's principle of 
optimality (ref. 10) it can be said that "any portion of an optimal trajectory is 
also optimal." However, the converse of piecing together portions of optimal tra- 
jectories will result in an optimal trajectory only under rather special conditions. 
But this piecing is precisely what we intend to do. It will be shown, from the 
numerical results, that piecing together a straight-line, optimal, high-altitude tra- 
jectory with a turning low-altitude trajectory will be an extremely close approxima- 
tion of an optimal flightpath. 

The first approximation is that the high-altitude portion of a long-distance 
optimal flightpath is a straight line in the horizontal plane. In reference 6 it was 
shown that for optimal horizontal flightpaths that involve initial and final turns, 
the flightpath does not contain a mathematically straight-line portion. However, even 
for flightpaths as short as 20 n. mi. there is a center portion of the flightpath 
that is straight for all practical purposes. Then, for the much longer flightpaths 
considered here, the assumption of a straight-line horizontal projection of the center 
portion of the trajectory appears reasonable. If, in addition, the initial portion 
of the turning descent is also almost a straight line and the speeds and controls at 
the junction match, so that there is no discontinuity, then one can be reasonably 
assured that the joining of the path segments will result in a close approximation of 
a fuel-optimal flightpath. 

In order to develop the performance index, Fj, it is convenient to refer to 
figure 1 where the flightpath has been drawn to include a turning climbout segment 
from tj to 10,000-ft altitude (t 2 ), a straight-line climbing segment (H-CLIMB^* from 
t 2 to cruise altitude, a straight-line cruise segment (RCRUISE) » a straight-line 
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descending segment ( r DESCENT) from cruise altitude to 10,000-ft altitude ( 1 3 ) » and a 
turning descent from 10,000-ft altitude to the final approach fix ti*. For a given 
stage length R, the optimal cruise altitude and speed change very little with minor 
changes in the length of the straight flightpath RcRUISE + RCLIMB + R DESCENT that 
may occur when considering alternative low-altitude climbouts and descents. Any 
change is primarily reflected in the change of the cruise distance, while the high- 
altitude climb and descent strategies are not affected. This means that the optimal 
fuel consumed during the straight-line climbout and descent segments is independent 
of the rest of the flightpath and can therefore be represented by a constant and the 
fuel consumed during the cruise segment is directly proportional to the range traveled 
at cruise altitude; then the total fuel consumed during the flight can be expressed 
by the equation 


F 


T 



F dt + cr cruise + 



F dt + Kj 


(7) 


where F is fuel burn rate, C is the fuel per nautical mile expended at optimal 
cruise speed and altitude (ref. 1 ), K 1 is the fuel consumed during the optimal 
straight-line climb and descent segments of flight, and ^CRUISE distance 

traveled at cruise. From figure 1 it is seen that 


r cruise r ~ r climb Descent x i X2 


where R (the distance between the climbout fix and final-approach fix), RcLIMB* and 
^DESCENT are constants. (Referring to fig. 1, we realize that y l and y 2 have 
negligible effect on the length of the long-distance straight-line flightpath above 
an altitude of 10,000 ft.) We define x x and x 2 to be: 

t 2 

v cos ^ dt ; x 2 
l 

so that 

J t 2 4 

(F - Cv cos ^)dt 4- I (F - Cv cos i|/)dt + K 1 + K 2 

Z 1 l 3 

where 



r 


v cos dt 


K 2 c(r r climb ^escent^ 

Therefore, provided that the overall flightpath is long enough to have a cruise 
segment, the performance function applicable to the low-altitude climb /descent por- 
tions of the complete flightpath has two parts and equals the fuel expended in the 
ascent (descent) minus the fuel saved by shortening the cruising portion of the high- 
altitude flightpath: 
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J 


(F - Cv cos ifOdt 


( 8 ) 


I 


,t? (tu) 


tl (t 3 ) 

The upper and lower limits of integration are those 
t 2 ), and descent (t 3 and t 4 ), respectively, and the 
free. 


appropriate for ascent (t x and 
upper limit of integration is 


The initial and final states for descent are: 


h 0 = 10,000 ft 

*0 = o 

v(0) < 250 knots IAS 
h f = 2000 ft 
= 180 knots 

and t|/f is specified between 0 and tt. The values are reversed for ascent. The 
values for (x 1 ,y 1 ) and (x 2 ,y 2 ) of figure 1 are not specified. The terminal altitude 
hf and speed Vf are chosen to permit a final straight-in landing approach along a 
specified glide slope. The initial speed v(0) is either determined by the speed 
constraint or, if the speed constraint is not violated, by the speed dictated when 
flying the solution of the optimization problem. Numerical results show that the 
terminal conditions at 10,000-ft altitude are the proper initial and final conditions 
for the straight en route portion of the flightpath. 


In order to solve the minimum-fuel climbout and descent problems it is necessary 
to model the drag, thrust, and fuel-flow rate for a specific aircraft. Choosing the 
Boeing 727-100 as an example, the drag, D, is modeled by 

D = T> 1 + D 2 u 2 (9a) 

where the drag for nonturning flight, D lf is given by 

D 1 = k 3 (l + k 2 h) + k 3 (1 + k 4 h)v 2 + k 5 (l + k 6 h)/v 2 (9b) 

and the bank-angle-induced drag, D 2 > is given by 

D 2 = ki(l - k 2 h) /2 + k 5 (l + k 6 h)/v 2 (9c) 

This approximation for drag results in a minimum drag of 9000 lb and includes 
sufficient flap deployment at the lower speeds such that an 8° angle of attack is not 
exceeded. In general, takeoff and landing weights differ from the nominal 
W = 150,000 lb. In this case, k l is multiplied by (W/150,000) and k 5 by 
(W/150,000) 2 to provide a reasonable approximation of the drag. 

The thrust value T max was chosen to be 23,000 lb and was set equal to 

zero. Both the maximum and minimum or idle thrust values are altitude- and speed- 
dependent. The idle thrust is small enough to be negligible. To simplify this 
analysis, the maximum allowable thrust has been chosen as the smallest maximum value 
allowable over the altitude range of 0 to 10,000 ft. 
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The fuel flow rate for the 727-100 is modeled by 

F = C 0 + C X T + C 2 T 2 (10) 

where 

Cj - C M+1 (1 + c, j+! h) + c tJ+s (l + c M+t h)vi j-o, 1 , 2 (11) 

and the numerical values of the k ? s and c f s are provided in appendix A. 

The optimal cruise fuel used per nautical mile, C, was set equal to 
17.5 lb/n. mi., which was obtained from reference 1. This particular C is for the 
example of a 200-n. mi. flight with a takeoff weight of 150,000 lb. This number will 
change for other weights and distances. We are interested here only in the character 
of the optimal turning climb and descent, and therefore do not change C or W in our 
computations. However, if optimization of a specific path was required, one would 
first optimize the straight-line path between airports and obtain C as well as 
takeoff and landing weights. This information would be accurate enough to compute 
the turning terminal-area paths without resorting to iteration. 


THE NECESSARY CONDITIONS 


In this particular optimal control, which requires the capturing of a specific 
heading, the state variables x and y in equations (1) and (2) can be ignored (hence, 
we do not have A x and Ay in the Hamiltonian); they are required, of course, for the 
flight trajectories in trie horizontal plane. 

For our analysis, we make the following assumptions, which are true for the 
numerical values of this problem. 

Assumption 1 . For the applicable range of velocities, available thrust exceeds 
the aircraft drag. This can be seen from figure 3 of appendix A where 
0 < D(v,h) < T m^ x . 

Assumption 2 . The fuel flow rate at T m i n =0 is positive, which eliminates 
the consideration of gliding flightpaths with shut-off engines. Also, 

T min > ~ C 1 (v,h)/[2C 2 (v,h)] for all v and h. This is significant for the calculation 
of the extremal thrust. 

The Hamiltonian for this problem is 

H = F + X T f + TiS 

= C 0 + C X T + C 2 T 2 - Cv cos + A^vy “ ^g u / v (12) 

+ X v g(T - D, - D 2 u 2 - W y )/W + n[v - v x (l + qh)] 

where n > 0, n(v - vj(l - qh) = 0 (ref. 11), and the convention of normalizing the 
costate variables by setting the costate variable, A 0 , which multiplies the perfor- 
mance function, to one is applied. The the remaining costate variables are given by 
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( 13 ) 


~ ^ 1 h'^ ~ + ^ 2 h u ^ ^ + r|V x c l 

\ = " H i|< = Cv sin * (14) 

*V = _H V = " C ° V " C lv T _ C 2 v T 2 " C COS * - V " V U/v2 

+ x v g(D lv + D 2v u 2 )/W - n (15) 

where = 9H/3h, and C 0 , = 9C 0 /Dh, etc. Since tf is not prescribed and H t = 0, 
then , 

H E 0, for all t e [0, t f ] (16) 


It is noted that in our approach the cruise cost C is part of the Hamiltonian, 
due to the performance index. The Hamiltonian would be identical if the fuel consumed 
had been the only variable used as performance index and if the values A x = C and 
Ay = 0 had been chosen. This observation ties the present work to reference 1, where 
the cruise cost A is the only adjoint variable and A* = C is the optimal cruise 
cost for a particular flightpath. 


Finding the argument of H, namely T, that minimizes H while using the per- 


where 


yields the extremal 

4 

thrust, T* = arg min H 
T 



T if 

max 

T > T 

max 


T* = « 

t if 

T . < t < T 

mm max 

(17) 


T . if 

mm 

T < T . 
mm 



T = ~(C 1 + 

A v g/W)/(2C 2 ) 

(18) 


and C 1 and C 2 are functions of h and v. Since minimization of H yields T* 
uniquely, it can be shown (see ref. 11) that T* and A v are continuous at junction 
times between the velocity-constrained and the unconstrained arcs. Thus thrust is 
seen to be a continuous function of A v and t. Since C 2 is small for all h and v, 
the range of A v for intermediate thrust is narrow and, as we shall see in the 
Results section of this paper, intermediate thrust will not occur for the optimal 
turning descents and climbouts. 

We observe that equations (4) and (14) imply that 

u = 0 if A^ = 0 on an interval (19) 

This is true irrespective of the minimization of H with respect to u. From the 
latter we obtain, by using H u = 0 and H uu > 0, 
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V 


if 


v 


and 


u 


* 


< u 


in 


X < 0 
v 


u sgn V if v > u and X < 0 
m ° 1 1 m v 


where 


v = -(X^W)/(2 X v D 2 v) 


( 20 ) 


( 21 ) 


If X v ^ 0 and X^ does not vanish on an interval (denoted by X^ = 0) , the mini- 
mization of H gives 


u* = u m sgn X^ ; X y > 0 and X^ ^ 0 (22) 

We note that X v cannot be positive while X^ vanishes on an interval because then 
minimization of H implies u* - ±u m , which is incompatible with equation (19). 
However, if X^ vanishes on an interval and X v crosses from negative to positive 
values, say at t = t 2 * then u switches from u(t) = 0, t < t 2 to u(t£) = ±u m , 
dependent on X^, which is no longer zero. This is a transition from a straight-line 
flightpath to a curved one. Such discontinuous bank angle would have to occur at 
zero thrust, thus it might occur (but actually did not occur) after initiation of the 
optimal descent. 


Lastly we evaluate the flightpath angle y. Since we have chosen to work with 
small angles (sin y = y, cos y = 1), y appears linearly in the equations. Therefore, 
y will either be at its limits based on the minimization of H, from y* = arg min H 


Y 

max 

if 

h y 

< 0 

Y . 
mm 

if 

h y 

> 0 


(23) 


where 


H Y = X h v - X vS 


or 


= X^v - X^g = 0 => X^ = X^g/v on a subinterval (24) 

in which case y will form a singular arc. In addition, when the speed constraint 
is reached, y is used to control speed. These latter two cases, where y will form 
a singular arc and where velocity is on the speed constraint, are covered below in 
more detail. 


Flightpath Guidance Along a Singular-y Arc 

When flying along a singular-y arc, equation (24) implies that all time deriva- 
tives of Hy vanish. 

(H y ) = A h v + X h v - A v g = 0 (25) 
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Solution of the simultaneous equations for H = 0 and Hy = 0 permit the calculation 
of the adjoint variables X v and A^. For simplicity, A v and A^ are calculated under 
the assumption (to be verified) that T = T max = constant for optimal climbout and 
T = Tmin = constant for optimal descent. There are three cases for u: u *= 0, 
u = u m (maximum), and -u m < u < u m (intermediate). We then examine the equations 
numerically to check whether they are consistent with the assumption (e.g., A v must 
be in the appropriate range for the assumed T, v must be in the speed range of 
interest, and the A f s must be real numbers). The singular arc flightpath profiles, 
are calculated for those cases where the singular arc is found possible. The 
singular arc can be found as follows. By using (24) in the equation for H = 0 (12) 
the dependence of (12) on y is eliminated. The additional assumption is made that 
the solution is not bounded by the speed constraint, ti = 0. In addition, by substi- 
tuting A^ from (13), A v from (15), v from (5), and A^ from (24) into (25), and 
by placing H = 0 and Hy = 0, we have three equations with three unknowns u, A v , 
and A^ . 

2H/v =► A (A + B u 2 ) - A . C u 
v v x x ip x 

h y / s - h/v =» x v (e x + p x u 2 ) + y: x u 


u l <u m J 


+ D x = 0 


-f G x = 0 


u = H A./A 

v I ; V 


(26a) 


with the following coefficients 

Aj . = 2g(T - D x ) / (vW) 

B x = -2gD 2 /(vW) 

C x = 2g/v 2 

D = 2F/v - 2C cos 
x r 

E x = °ih v/w ■ e D iv /w 
P x = D 2h v/W - D 2v g/W 

G x = ~ F/v - F h v/ S + F v 
H x = -W/(2D 2 v) 

Solving equations (26) for A v and X^ with |u| < 1 % 


(26b) 


L * [D (P H + C ) - G (B H 

v XXX X X XX 


- c )]/[E (B H - C ) - 

X 7 1 1 X X X X 


A (P H + C )] 

XXX X 


(27) 


= 


X /(G A - D E )/{H [D (P H + C ) - 

C X X X X X X X X X 


G (B H - C ) ] } 

XXX X 


(28) 
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In the special case u = 0, from (19), 


A , = 0 


(29) 


and (26) becomes 


A A 

+ D =0 

(30a) 

V X 

X 

A E 

+ G v = 0 

(30b) 

V X 

X 

1 A 

- D E =0 

(31) 

X X 

X X 

, = o, 

. The variables G x , A x , D x , and E x 

are 


which to be consistent requires 


This is consistent with (28) for 
functions of altitude and speed so that we can expect a single altitude-versus-speed 
profile. If we attempt to obtain an explicit expression for h versus v by expanding 
(31) we encounter an extremely long polynomial in the 7th power of v and 2nd power 
of h, which is not very illuminating. We therefore shall be satisfied with a numeri- 
cal solution of (31). In this case, the altitude-speed profile defines the state 
space and we do not need to solve Hy = 0 to determine y. Instead, by replacing h 
with Ah/ At in (3) and v with Av/At in (5) and by dividing the equations and 
solving for y, y can be determined directly from the altitude-speed profile. 


y = [(T - D 1 )/W][(Ah/Av)/(Ah/Av + v/g) ] 
Solving equation (26) for X v and A^ with |u| = u m 


(32) 


= -(D + G ) / [ (A + B u 2 ) + (E + P u 2 )] 

x x x xm x xm 


(33) 


and 


X ,= [X (A + B u 2 ) + D ]/(C u ) 
<f> v x x m x x m 


(34) 


This singular case is of little interest for optimum ascents and descents, since 
u = %ax was almost never reached on the singular-y arc. It should be noted that 
while A v and A^ can be computed directly from the above equations all along the 
singular-y arc, it is computationally simpler to use the above equations for computa- 
tion of initial values only and to integrate (15) to obtain A v and (14) to obtain 


The six cases, ascents and descents, each with u = 0, u = u m , and -u m < u < u m , 
must be checked to determine whether singular arcs are possible. Since the equations 
are altitude dependent, these cases are checked numerically at h = 10,000 ft, which 
is useful for forward integration in descent and backward integration in climb. From 
the numerical solution of (31), A v is obtained from (30a) or (30b). This A v , used 
in (18), is consistent with T = T^x for climb and T = 0 for descent. For descent 
the solution of (31) gives a speed below 250 knots IAS, which means that the speed 
limit will be inactive as required. During the ascent the speed profile along the 
singular arc crosses above the speed limit. Thus, strictly speaking, ascent on the 
singular-y arc is not permitted. However, this case is investigated, since the speed 
limit in the terminal area may be relaxed by ATC. For u = u max , A v is consistent 


11 



with the assumption for the thrust. Thus singular arc profiles are possible for all 
u. The above conclusions are identical for all altitudes. Being assured that the 
singular-Y case is indeed a candidate for an optimal flightpath, the singular arc 
flightpath for intermediate u is determined from the condition (Hy) = 0: 

0 = (Hy) = 3(Hy)/3t = 3/3t{-F h v + A^DjV + g(T - D)/v - gD v ]g/W 
+ g(F v - C cos ip + A^gu/v 2 ) } 

= -F h^ " F h V + VV + g(T ‘ D)/v " SD v ]g/W + A v [D h v + D h v - g(T - D)v/v 2 1 (35) 

- g(° u u + D h k + D v v)/v - gD v ]g/W 

+ g(F v + C sin ipty + A^gu/v 2 - 2A^guv/v 3 + A^gu/v 2 ) 

For definition of the fuel flow and drag derivatives see appendix B. To solve 
(35) for y is straightforward but tediousand will only be.sketched out. First, all 
terms with different time derivatives — v, h, i p, A v , X^, u, F v , F^, D v , and are 
combined. The last five time derivatives are eliminated with expressions (B13), (B3) , 
(B5), (B9) , and (B 11) from appendix B. This results in the following equation for 
(Hy) 

<V - V + B H h + C H * + D H i v + (36) 

where the coefficients have the following values: 

A h = -F h + gD h A v /W - g 2 A v (T - D) / (Wv 2 ) - 2g 2 X |J u/v 3 - g 2 A v D v /(Wv) 

- Q H (l/v + D 2 v /D 2 ) - vF vh + 2gk 3 k 4 v 2 X v /W - 2gk 5 k 6 (l + u 2 )X v /(Wv 2 ) 

- 2g 2 k 3 ( 1 + k 4 h)X v /W - 6X v g 2 k 5 (l + k 6 h)(l + u 2 ) / (v**W) 

B H = -g 2x v D h /( Wv ) - Q H D2 h /D2 + gF vh " 2 g 2k 3 k 4 vX v /W + 2 g 2k 5 k 6 (1 + u 2 )X v /(Wv 3 
C R = gC sin 4) 

D r = g(D h v + g(T - D)/v - gD v )/W - Q h /X v 
E r = g 2 u/v 2 + Q h /X^ 

Q h = t g 2 A^/v 2 - g 2 A v D u /(Wv) + (k x k 2 v + 2k 5 k 6 /v)gA v u/W 

+ 4g 2 k 5 (l + k 6 h)uX v /(Wv 3 ) ]u 

None of these coefficients involve y . When the equations for the time deriva- 
tives are inserted in (36), y appears linearly and the solution for y is 
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Y = [AjjgOT - D)/W - C R gu/v + D r (-F v + c cos f - gX^u/v 2 + 8 X V D V /W) 
+ E r (-Cv sin ] / [A R g - B R v + D R X h ] 


( 38 ) 


For (38) to be a candidate for an optimal solution the generalized Legendre-Clebsch 
condition must hold 

<-H y ) y > 0 * <-H y ) y - - B h v + D H X h > 0 (39) 

The condition (-Hy) =0 is not possible since (39) is also the denominator for (38) , 
which would imply y 00 . 

Again we emphasize that the above development applies to the case of constant 
thrust only since 3T/3t = 0 was assumed. 


Flightpath Guidance Along the Speed Constraint 

Since the climb speed for the singular y was above the speed constraint we must 
now develop the equations for y when flying at the speed constraint, n ^ 0. Equa- 
tion (6) is solved for v, which is then differentiated with respect to time and h 
is replaced with (3). Setting the resulting expression for v equal to (5), and 
solving for y results in the following expression 

y = [(T - D)/W]{g/[g + v 2 q(l + qh) ] } (40) 


where the second factor varies between 0.92 and 0.91 over the altitude range of 
10,000 ft. Since S = 0 on the speed constraint, equation (6) can be solved for v, 
and the derivative taken in order to express v and h on the constraint 

v = Vj(l + qh) (Ala) 

v = v^.qh (Alb) 

but 

h = v^.(l 4- qh)y (Ale) 

which implies 

v = qv 2 (l + qh)Y (42) 

Equations (13) and (15) can now be solved for n* In addition, we must satisfy 
Hy = 0, since y was not determined by this necessary condition. The solution pro- 
ceeds as follows 


H = 0 = A ■ v - A g A, = A g/v 
y h v & h v & 


differentiating equation (A3) results in 


= g(vX v - X v v)/v 


(43) 


(44) 
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Combining equation (44) with equations (13) and (15), we can obtain an expression for 
n which is not dependent on X^ or X v 

I [ C o h + C lh T + C 2 ^T 2 - X v gD h /W - g(vA + X v v)/v 2 ]/(v ]; q + g/v) >0 on the 

41 1 1 constraint 

0 off the 

constraint 

(45a) 

where 

A n " C 0v + C lv T + C 2 v T 2 - C cos * + X h y + X^gu/v 2 - X v gD v /W (45b) 

In (A5) v and v are given by (41a) and (42). 

To complete the statement of the necessary conditions, the character of the 

velocity set must be studied. The above results may be suspect, since the velocity 
set for this problem (v, h, F) is not convex everywhere. However, it is shown in 
appendix C that the nonconvexity of the velocity set may be ignored for this minimum- 
fuel problem. 


COMPUTATION OF EXTREMALS 


The extremals are computed by numerical integration of the state and costate 
equations ( 1 ) — (5 ) and. ( 13)— ( 15) with the control T given by (17)-(18), u given by 
(19)-(22), and y given by (23) or (38) on the singular-y arc or by (40) when the 
trajectory is on the speed constraint. On the singular-y arc, the costate variables 
X v and A^ are obtained by integrating (13)— (15) instead of finding them directly from 
(27) and (28); we find A^ directly from (24), whereas on the speed constraint 
Aft is computed directly from (43). In this relatively simple case of a heading 
capture with specified altitudes and speeds, the costate variables can be treated as 
parameters in such a manner that families of extremals will sweep out all desired end 
conditions. To achieve this, forward time integration must be used for descent and 
reverse time integration for ascent. For convenience of integration the (x,y) coor- 
dinates at the 10,000-ft altitude were chosen (0,0). This technique permits the path 
integration to be started on a constrained or a singular portion of the path, which 
limits the freedom of choice for the adjoint variables and makes it easier to pick 
them from the limited set. When starting on a singular portion of the path, the opti- 
mal path leaves the singular arc as required by the end conditions, provided the 
solution is not forced off by failing to meet the speed constraint or the generalized 
Legendre-Clebsch condition (eq. (39)). When starting on a speed-constrained portion 
of the path, the optimal solution will leave the constraint boundary as required by 
the end conditions, provided the solution is not forced off by failing to meet n > 0. 
If in the singular-y-arc descent, the speed constraint is violated before reaching 
the reference altitude of 2000 ft, the guidance switches to the speed-constrained arc 
at this point. The singular-y case and the speed-limited case are discussed 
separately. 


14 



The Singular-y Case 


As previously discussed, A v and A^ can be obtained from (27) and (28) and A^ 
from (24). With the thrust T = T max or T = T in i n and h given (h = 10,000 ft), 
the only remaining variable of choice is v. In addition, A^ must be close to but 
not equal to zero for turning to occur in the lower-altitude portion of the flight- 
path with almost a straight-line flight in the higher-altitude portion. Therefore, 
we choose v iteratively to solve (31), which is the condition for A^ = 0 in (28). 
Small variation from the speed thus found will permit us to generate extremals which 
cover the range of heading 0 < ip < 180°. Changes in v to cover the range of head- 
ings were so slight (in the order of 10” 19 knots) that double-precision calculations 
were required. For all extremals, the generalized Legendre-Clebsch condition (39) 
remained valid over the range of altitudes considered. For practical purposes, we 
determined v from (31), A v from (30a) or (30b), and selected A^ over a small 
range to cover all turns (0 - 180°). 

Based on the constraint for flightpath angle, y max = in descent anc * 

Ymin = °° in an initial descent in the optimal climb trajectories or a final 

climb in the descent trajectories are not allowed, even though they may have lower 
cost. If in the singular-y-arc descent the speed constraint is violated before reach- 
ing the reference altitude of 2000 ft, the guidance must switch to the speed- 
constrained arc. In most cases, getting off the singular-y-arc or the speed- 
constrained arc at the lower reference altitude h = 2000 ft transfers the path to 
the zero degree glide slope angle and maintains it there while turning and reducing 
the speed to 180 knots. For optimal descent with large turn angles, the optimality 
condition (23) demands a step change of glide slope angle which brings the extremal to 
an incorrect final altitude. By getting off the singular-y arc at a higher altitude, 
the proper final altitude at the proper speed could be obtained after a few itera- 
tions. Whatever final heading is obtained, the resulting extremal provides one 
example of an optimal heading change climbout or descent maneuver. A second method 
to force a flat turn after leaving the singular arc consists of constraining the 
permissible y to a narrow range. 


The Speed-Limited Case 

If A v is treated as a parameter, A^ is determined via (43); Aw, can be 
obtained by solving equation (12) and replacing A^ with equation (4 j) , u with 
equation (21) for an unsaturated bank angle, and y with equation (40). 


A^ = ±2v/-[C 0 + C X T + C 2 T^ + Cv - A v g(T - J ) 1 ) /W] A y D 2 / (gW) 


(46) 


Again, the choice of one adjoint variable determines the other two once the thrust is 
chosen (T = 0 or T = T max ) . As before, initially u, and hence A^, must be 
extremely close to zero. Hence, from (46) A v must be close to 


A 


v 


A 


l|/=Q 


= -(C 0 + CiT + C 2 T 2 - Cv)W/[g(T - D 2 )] 


(47) 


Now from (18) and (20), A v must be in the range 

-WC^/g < A v < 0 for T = 0 


(48) 
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and from (18) 


X < - (C T + 2CLT )W/g for T = T (49) 

v 1 2 max & max 

The above considerations limit the possible values of A v that can be considered 
when starting on the speed constraint. In addition, n in (45) must be positive, so 
that an extremal on the constraint can be considered a candidate for an optimal 
f lightpath. 

We remember that C in (47) is the cruising efficiency: pounds of fuel per 

foot traveled at the cruising speed. The number depends in a minor sense on the 
stage length. Or, for C = 0 in the performance criterion, we would obtain the 
minimum fuel climbing or descending heading capture flightpath, which does not con- 
sider any continuation of the path. Therefore, the permissible range of A v shall be 
translated into the equivalent range for C. For descent, T = 0, the two limits in 
(48) for A v are used and solved in (47) for C. This gives 

(C 0 + C 1 D 1 )/v > -C > C 0 /v 

numerically 

0.00407 > -C > 0.000148 

The equation for n(0), (45), is a linear function of C and must be positive. 
Numerically, n(0) just passes through zero as C becomes greater than -0.00237. 

Hence, C must be in the range 

0.00237 > C > 0.000148 

Hence both C's that are considered (C = 0 and C = 0.00288 Ib/ft = 17.5 lb/n. mi.) 
are outside the allowable range. It is therefore concluded that descent on a 250-knot 
speed bound is not optimal. This conclusion will be confirmed by results for the 
singular-y arc where it has already been shown that the speed of descent on the 
singular arc is below the speed constraint of 250 knots. Similar considerations for 
the climbout show that extremals calculated for both C's are candidates for optimal 
flightpaths. Again, small variations from the value given in (47) will cover all 
possible heading changes (0 - 180°) and as in the singular-y case the sign ambiguity 
for A^ (46) provides for right and left turns. 


RESULTS 


We will first discuss the special case for straight-line ascent and descent 
(A^ = 0) to tie our data with that of (1). This will be followed with results for 
optimal turning ascents and descents. 

Figure 2 is a plot of the altitude-versus-speed profiles for straight-line climb 
and descent. The thin lines are a reproduction of a figure from (1) for a 200-n. mi. 
altitude-speed profile, which we will use for verification of our results. The 
singular-y-arc ascent profile, denoted by ' 1 ? in figure 2, is above the FAA terminal- 
area speed constraint. The altitude-versus-speed profile connects to the profile for 
the complete 200-n. mi. path (point A) at a speed somewhat less than the speed at 
which the ascent is resumed from the 10,000-ft altitude (point B). This is a result 
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1 - SINGULAR ASCENT 

2 - 250 knot INDICATED AIRSPEED ASCENT 
3-MINIMUM FUEL CLIMB (NO CRUISE) 



TRUE AIRSPEED, knots 


(a) Altitude-speed profile. 



FLIGHT PATH ANGLE y, deg 

(b) Flightpath angle command. 

Figure 2.- Straight-line optimal flightpaths. 
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of a conservatively chosen T max of 23,000 lb. By setting T max at 25,000 lb the 
10,000-ft intercept shifts to 338.5 knots, indicating that the following climb is 
the continuation of the singular-y arc. It is thought that this is closer to the 
effective T max used in (1). As discussed before, the FAA often permits a climbout 
faster than 250 knots IAS. In this case, the optimal (singular-y-arc) climbout saves 
36 lb of fuel for the total path compared with the speed-limited climbout, denoted by 
f 2’ in figure 2, which is followed by an acceleration to the singular-y-arc speed at 
the' 10,000-ft altitude. Figure 2 also shows (dash-dot line f 3 T ) the singular-y-arc 
climb for C = 0, which would be the optimum climb profile if the cruise did not have 
to be considered. Since an acceleration to the speed at point B (fig. 2) is required, 
the total cost is higher than for the singular-y-arc climb with the proper C, but 
less than the cost for the speed-limited climb. 

The singular-y-arc optimal descent speed, denoted by ’A’ in figure 2(a), is 
below the speed limit of 250 knots IAS, denoted by 1 2 1 in figure 2(a) and turns more 
sharply toward the speed constraint than the altitude-speed profile given in (1). 
Therefore, it approaches the speed constraint at a higher altitude than the example 
of (1). The altitude-speed profiles are very sensitive with respect to small differ- 
ences in the engine-idle fuel flow model. For instance, when the modeled engine-idle 
flow is reduced to one-third of the present value, the singular-y descent altitude- 
speed profile will always remain below the 250-knot indicated airspeed constraint. 

Once the model is fixed, the resulting fuel consumption is not very sensitive with 
respect to the altitude-speed profile. For instance, flying the mixed descent or a 
pure singular-y-arc descent (dashed-line extension of ’A T ) makes less than a pound 
difference in fuel for the straight-line or even turning descent. Figure 2(b) also 
shows that for both climb and descent, the y's remain almost constant, which makes 
flying these flightpaths relatively easy. With our results for straight-line ascent 
and descent agreeing well with those of reference 1, we shall first look at the data 
for turning-ascending flight and then at those for turning descents. 

Solutions to the optimization problem that include ascending turns result in 
altitude-speed profiles that are almost identical to the straight-line optimization 
problem (see fig. 3(a)), with most of the turn occurring in level flight before ascent 
from the minimum maneuvering altitude (2000 ft in this example; see fig. 3(b)). To 
further illustrate the above statement, table 1 summarizes the numerical values for 
the key flightpath variables as a function of the total heading change for both 
singular-arc and speed-constrained ascents. We notice from table 1 that the speed at 
which ascent begins (column 3) depends on the type of ascent and not on the required 
total heading change (column 2). A similar statement is largely true for the initial 


TABLE 1.- NUMERICAL VALUES OF KEY FLIGHTPATH VARIABLES IMMEDIATELY 
AFTER SWITCH-OVER FROM LEVEL TURN TO CLIMB 


Case 

(i) 

^total ’ 
deg 

(2) 

V, 

knots 

(3) 

Y at 

h = 2000 ft, 
deg (4) 

to go, deg 
(5) 

<f> . deg 
(6) 

Type of 
ascent 

I 

0 

301.7 

4.356 

0 

0 

Singular arc 

II 

90 

301.7 

4.355 

1.6 

1.7 

Singular arc 

III 

180 

301.7 

4.353 

5.2 

5.2 

Singular arc 

IV 

0 

258.1 

4.79 

0 

0 

Speed constrained 

V 

90 

258.1 

4.76 

9.2 

7.2 

Speed constrained 

VI 

180 

258.1 

4.65 

30 

17.8 

Speed constrained 
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n.mi 



x, n.mi. 


(b) Horizontal projection. 

Figure 3.- Optimum turning climbout on the singular-y arc. 
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K ANGLE <t>, deg 


BANK ANGLE LIMIT 


BANK ANGLE LIMIT 



Av!/ = HEADING CHANGE TO GO, deg A^ = HEADING CHANGE TO GO, deg 

Figure 4.- Bank angle versus heading-change-to-go in optimal climbout. 

= maximum heading change. 


1 - OPTIMUM CLIMBOUT 
FOR MINIMUM FUEL 
OVERALL TRAJECTORY, C ¥= 0 



x, n. mi. x, n. mi. 


(a) Horizontal projection. (b) Bank-angle histories. 

Figure 5.- Comparison of two types of optimal climbouts. 
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, deg 


flightpath angles, which vary only minimally (column 4). The bank angles experienced 
during the entire turning singular arc and speed-constrained climbouts are shown as 
functions of the heading change in figures 4(a) and 4(b), respectively. We notice 
that the bank angles gradually go to zero as the desired heading (Ai|/ = 0) is 
approached. As table 1 shows for our example (case III), the bank angle at the point 
of switching from level flight to climb is already as small as 5.2°, indicating that 
the ascent is almost a straight line. Also, for a turn of more than 90°, the early 
part of the turn is flown nearly at or at the bank-angle saturation. Intuitively 
this makes sense, because while flying away from the target we must turn as fast 
as possible. (This is true in spite of the opposite effect that maximum bank angle 
has on acceleration, which also must be maximized.) The penalty of observing the 
speed constraint is about 36 lb of fuel, independent of the amount of turn. 

The altitude-speed profile for the singular-y-arc climb (fig. 3(a)), can be well 
approximated by a constant IAS speed climb. If we set the speed constraint to a 
value that corresponds to the singular-arc speed at 10,000-ft altitude — 331.4 knots — 
the resulting performance *is basically identical with the singular-y-arc performance. 

In figure 5 we compare the minimum fuel climb, which considers the fuel saved 
by flying a shorter cruise distance with the equivalent minimum fuel climb (C = 0) 
for a heading change of 68°. The latter indeed takes 7 lb less fuel, but it climbs 
in an x-distance that is 2.93 n. mi. shorter, thus requiring an additional 52 lb of 
fuel for the longer cruise. The bank-angle histories, shown in figure 5(b), are 
quite different for the two cases. 

For the optimum descent only an extremal starting with a singular-y-arc descent 
is a candidate for optimality (fig. 6(a)). We note here that this turning descent 
has actually three phases: singular arc, constant indicated airspeed arc, and level 

deceleration. Again, most of the turn occurs in level flight (fig. 6(b)) and the 
altitude-speed profiles are identical within the power of resolution of the graph 
(fig. 6(a)). Figure 7 shows the bank angles versus heading. The bank angles essen- 
tially follow the same curve independent of the total heading change. This allows a 
simple approximation of the optimal descent in an onboard algorithm. Note that for 
this example the beginning of the turn is executed with gradually increasing bank 
angle until, for heading changes greater than 25°, the bank angle is at its maximum. 
For turns greater than 90° this is required to turn as fast as possible with as large 
a deceleration as possible to the final speed and heading. The performance index 
increases with the turn angle, since only the x-distance flown reduces the en route 
portion of the flightpath. However, the fuel for the final descent decreases slightly 
for larger turn angles, since the increased bank angle increases drag, which brings 
the aircraft more quickly to the desired final speed. 

Table 2 summarizes the numerical values for the key flightpath variables as 
functions of the total heading change, The entries in table 2 are similar to 

those in table 1, except for the addition of column 3, which indicates the altitude 
at which the system leaves the singular-y arc. The situation for descent is somewhat 
more complex than for the ascent; this will be explained by means of the six represen- 
tative cases in table 2. For small turn angles (see cases I and II, table 2), the 
singular-y descent to 2000 ft is followed by a horizontal decelerating turn. 

As stated earlier, for large angles of turn the optimal path required a change 
in y from ymax (level flight) to Y m in P art wa Y in the fina l turn, after having 
left the singular-y arc (see cases III and IV, table 2). The changes in altitudes 
for y-switching that are required would make an onboard algorithm quite complex. 
Therefore, the following test was made. Upon leaving the singular-y arc we forced 
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ALTITUDE, ft 






BANK ANGLE LIMIT 


0 25 50 75 100 125 150 175 

A\p = HEADING CHANGE FROM INITIAL HEADING, deg 

Figure 7.- Bank angles for singular-y-arc descent. 
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TABLE 2.- NUMERICAL VALUES OF KEY FLIGHTPATH VARIABLES IMMEDIATELY BEFORE 

SWITCH-OVER TO LEVEL TURN 


Case 

(i) 

^T> 

deg 

(2) 

h switch 
to level 
turn 
(3) 

V, 

knots 

(4) 

YSA* 

deg 

(5) 

during 
descent , 
deg (6) 

deg 

(7) 

Comments 

I 



258.6 

-3.35 

0.0 

0.0 

Ymax = 0 » Ymin = 4 

II 

■ 

i 

258.6 

-3.35 

0.0 

.02 

System wants to stay at 

hi 

mm 

2025 

tm 

-3.28 

3.6 

3.6 

Ymax “ 0 , Y m i n = -4 

IV 

a 

2591 

m 

-3.28 

4.6 

4.6 

System switches to Ymin 
after initial level turn 

V 

mm 

— 

258.6 

-3.35 

1.3 

2.3 

Ymax = 0 * Y m i n = -0.001 

VI 

IB 

m 

258. 1 

-3.78 

20.0 

30.0 

System forced to level turn 


the aircraft to stay at 2000 ft by setting very narrow limits for maximum and minimum 
y (see cases V and VI of table 2), We find that there is a negligible difference 
between performances for the optimal case with the y-switching and the corresponding 
constant-altitude decelerating turn. Beyond = 170° a bank angle of 30° (u xnax ) 
is reached on the singular arc, which switches to a positive Ysa f° r a 
decelerating turn. Since Ymax ” 0° violated, we must get off the singular-y arc 
as soon as u max is reached. However, for purposes of a possible onboard algorithm, 
we studied a 180° turn singular-y descent such that when u m ^ x is reached before the 
2000-ft altitude, we then continue nonoptimally at constant Ymin = -4° to 2000 ft, 
followed by y max = 0° such that Vf± na i = 180 knots is reached at the completion 
of the turn. Again, there is only a 3.2 lb of fuel increase for this case compared 
to the optimal case. 

To summarize, for a simple descent, which involves a final level turn, the 
altitude-speed profile deviates only slightly from the straight-line flight even for 
the largest heading change (column 4 of table 2). Similar to the ascent, most of the 
turn is accomplished in level flight (^p - A\J> of columns 2 and 6 of table 2). 

The question of global optimality remains. We must find all extremal controls 
and choose one with the lowest cost. We have shown earlier that on the singular-y 
arc and on the speed limit we have single-parameter families of extremals. The varia- 
tions of those parameters were very small in order to cover the range of final head- 
ings. For larger changes we get descending or ascending turns that include a 360° 
loop. Since our object is to cover large distances in the x-direction, such large 
turns will not be optimal. Also, from the problem formulation, turns above 180° will 
not be optimal, since they could be replaced by smaller turns in the opposite direc- 
tion without increasing the cruise cost, and the cost of the turn always increases 
with the magnitude of the heading change. We are aware, however, that we have not 
studied extremals that do not include either a singular arc or a speed-limited sec- 
tion. Therefore, while we have not proved this, we have reason to suspect that our 
flightpaths are indeed globally optimal. 
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COMPARISON OF OPTIMAL PROCEDURES WITH VARIOUS 
SUBOPTIMAL PROCEDURES 


The question is often asked how much fuel is saved by using the optimal proce- 
dure. Such questions are difficult to answer, since we must ask: "What should the 

optimal procedure be compared with?" We shall compare reasonable suboptimal climbout 
procedures which more and more closely mimic the optimal. First let us assume that 
the optimal climbout altitude-speed profile for straight-line flight is known, and 
that we will fly at least the latter part at an approximation of the optimal altitude- 
speed profile. Such an approximation is a climb at maximum thrust at a flightpath 
angle, which is the average of the singular gamma angle. This assures us that we will 
meet the terminal conditions. In addition, for all climbout procedures starting at 
2000-ft altitude we shall execute an initial maximum bank-angle turn until the 
desired heading is reached and continue with straight-line flight. The only varia- 
tions in the climbout procedure considered are the initial climb angles until the 
singular-y speed is achieved. 

The results are shown in figure 8. Let us first consider an initial 3° climb. 

The additional fuel is between 10 and 28% over the optimum, which is 285 lb of addi- 
tional fuel for a 180° turn. If we climb at a slower rate of 500 ft /min, the extra 
fuel cost is about 4% or 36 lb for the 180° turn. The best approximation of an opti- 
mal turning climbout is an initial horizontal flight . Here the additional cost is 
primarily due to the incorrect bank-angle schedule. The smallest additional cost 
beyond the optimal is due to flying a constant y as an approximation of the 
slightly varying singular y. This is illustrated in figure 9, where we plot the 
additional fuel when we change the climb angle from the average singular y. Over a 
range of 0.04° the cost deviates from the optimal by only 0.1 lb, and it is fairly 
insensitive for errors of ±0.15° from the singular y. Considering these results, we 
can say that the most important fuel-saving feature of the optimal climbout is the 
initial horizontal turn and acceleration up to the singular-y speed. In fact, if 
an initial dive would have been permitted, more fuel could have been saved. 

We can make a similar investigation for the turning descents. Again, we descend 
on a flightpath angle, which is the average singular-arc descent angle. At a certain 
altitude we reach a speed of 250-knots indicated airspeed. At this point we select 
the varying flightpath angle which will hold this airspeed. When we reach the 
desired altitude, we go to the final flightpath angle, which is nominally 0 , and 
begin the turn at maximum bank angle in such a manner that the desired heading and 
final speed are reached simultaneously. 

It should be noticed for the following that the optimal flightpaths with which 
we are making the comparison are constrained to a level final turn. We stated 
earlier that for large turns, where the maximum y was constrained to 0° and the 
minimum to -4°, it was optimal to begin the level turn at a higher altitude than the 
final altitude and finish the turn at -4° to end up at the desired final altitude 
and speed. 

We compare the constrained optimal descent with various suboptimal ones, where 
the only change is the flightpath angle of the final turn (see fig. 10). For the 
descent that mimics the optimal descent most closely (final flightpath angle = 0 ), 
the suboptimal performance is within 0.2% for all heading changes. Again, a final 
descending turn will usually result in an increase of fuel used. For large turns this 
trend reverses, since this approximates better the less constrained final optimal 
turn. If we permit a final climb (and the initial undershoot in altitude, which is 


25 




0 30 60 90 120 150 180 

HEADING CHANGE, deg 


Figure 8.- Extra fuel expenditure for various suboptimal climbout strategies 



7, deg 

Figure 9.- Change in climbout fuel cost for various climbout angles. 


connected with it) for small angles of turn, the suboptimal descent performs slightly 
better than the optimal, which was constrained not to climb during the descent. Since 
descent is executed with idle thrust, the above changes in performance are not as sig- 
nificant in terms of fuel as the changes from the optimal for the turning climbout. 


CONCLUSIONS 


1. For onboard optimization of fuel-efficient turning descents and climbouts, 
it is sufficient to use an existing climb-cruise-descent fuel-optimization algorithm 
and superimpose initial- and final-level turns from and to the airport. 

2. Relaxation of the speed limit for climbout of the 727-100 can save 36 lb of 
fuel per climbout, independent of the required heading change. 

3. A constant-IAS climb after an initial-level turn using the optimal bank- 
angle schedule is a close approximation to the optimal turning climbout. 

4. An initial-level accelerating turn to the proper speed before climbing is the 
most fuel-saving feature in an optimal climbout. 

5. The optimal descent starts below the FAA terminal-area speed constraint of 
250 knots IAS, then reaches the speed constraint and finally terminates in a decel- 
erating level turn. 



Figure 10.- Optimal versus suboptimal descents. 
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APPENDIX A 


The Drag and Fuel Flow Equations 

Reference 1 includes subroutines that calculate drag thrust and fuel flow from 
comprehensive table lookup data over the complete operating range of the Boeing 
727-100 aircraft. For our method of calculating extremals, we need functional rela- 
tionships for the limited range of 150- to 350-knot speeds and 0- to 10,000-ft 
altitudes, which are derived here. 


Drag D(v,h) 


Over the range of speeds considered, flaps must be deployed for the lower speeds. 
However, we do not want to consider flaps as a control variable. Therefore, flaps 
are eliminated by putting a limit on the angle of attack 


a < a =8° 
“ max 


(Al) 


and by automatically deploying the minimum amount of flaps so that a = 8° is never 
exceeded. Although the drag curves are shown for all altitudes, in practice flaps 
are deployed at low altitudes only. This is in agreement with the optimal results. 
Since we have the constraint 


L = W/cos <f> 


(A2) 


the drag is also bank-angle dependent. 

D = D(v,h,<}>) 


(A3) 


It should be noted here that our flap schedule is different from that in refer- 
ence 1, which deploys flaps as a function of indicated airspeed, as shown in 
figure Al. Our flap program results show less drag at the lower speeds than the 
schedule of reference 1. However, neither method represents the real world com- 
pletely, since, at present, flaps are controlled manually in steps, based on speed 
cues . 


The calculation of the drag proceeds as follows. Given h and v from the 1962 
atmospheric model, we obtain density p and Mach number, M. Knowing the bank angle, 
we calculate the required lift coefficient to meet the constraint (A2) 

C L = W/(0.5 pSv 2 cos <p) (A4) 

where S = area of the wings = 1560 ft, W = weight of the aircraft = 150,000 lb, and 
p = density of air at the given altitude. Then, for a = 8° the subroutine CLIFT 
computes the available lift coefficient CL av (a,h,M, 6f ) | a=0 0 . Starting with zero 

flaps, the flaps are increased in 0.01° steps until 

C T > C T (A5) 
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The final flap value is the value used in the subroutine CDRAG which determines the 
drag coefficient Cj)(M,Cl , 6f ) . This permits us to calculate the drag 

D = 0.5 PSV 2 C D - D(h.v,«| ose „ jL ^ /co8t (A6) 

An example of the drag curves obtained by the above method, including the flaps 
required, is shown in figure A2. The bumpiness of the drag-versus-speed curves in 
the region where flaps are deployed is an artifact of linear interpolation of the 
drag polars which are given only for 0°, 2°, 5°, 15°, 25°, and 40° flap angles. Our 
model therefore should smooth out these bumps. 


Our task is to match the drag data with an analytic expression. In order to 
minimize the number of parameters required, we shall first develop an expression from 
the approximate theory of wings, and then determine parameters that fit the data. 


The lift coefficient of a clean wing is 

C L = C L ° 
a 

Adding flaps shifts the curve upward 

C i.' C L 0 +C L tl 
0 a 

Without flaps the lift-drag polar is 

S ■ C D 0 + ’ C D 0 + "V 

Assuming the same form when the drag polar is shifted due to flaps adds a linear 
term in a. 


(A7 ) 


(A8) 


(A9) 


c i - + y' c l + ' 


’( C L “ + C L '‘I 

\ a a / 


(A 10) 


From (A4) and (A8) we can solve for 

a = W/(0.5 pSv 2 cos <f>) - C T 1/C T 


(All) 


and using (All) to replace a in the drag equation (A10) results in an equation for 
drag of the form 


D = 0.5 pSv 2 = k 0 /cos <j> + k x v 2 + k 2 (l + tan 2 4>)/v 2 


(A12) 


where the k ? s are functions of p, S, W, n 1 , C^, anc * 

This is the form of the functional relationship for a constant altitude and 
constant flap angle. In actuality the flap changes with speed for constant altitude 
as shown in figure A2. The bank angle, which is one of the controls, appears twice 
in the above expression. The first approximation of the k’s in (A12), which must 
be functions of altitude, are found as follows. We write for 4) = 0 
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D = k x (1 + k 2 h) + k 3 ( 1 + k 4 h)v 2 + k 5 (l + k 6 h)/v 2 


(A13) 


Using the tabular data graphically represented in figure A2 , for h = 0, and for 
three different speeds of 150, 225, and 350 knots, we obtain three equations with 
three unknowns, k 1 , k 3 , and k 5 . With these solved, we chose h = 10,000 ft and 
solve for k 2 , k 4 , and k 6 . These six constants are a good first approximation to the 
drag equation. These constants are then used for the complete drag equation (A12) 

D = k 1 (l 4- k 2 h)/cos cj> 4- k 3 (l 4* k^h)v 2 4- k 5 (l 4- k 6 h)(l 4- tan 2 4>)/v 2 

which using cos $ = 1//1 4* tan 2 4> * 1/(1 4- 0.5 tan 2 tf>) is approximated as 

D = k^l + ^h) + k 3 (l + k 4 h)v + k 5 (l + k 6 h)/v 2 + [k x (1 4- k 2 h) /2 + k 5 (l 4- k 6 h)/v 2 ]tan 2 cf> 

(A14) 

For the purpose of an analytic approximation of the drag function, we generate 
180 samples for all combinations of the following parameters: 

4> = 0°, 10°, 20°, 30° 

h = 0, 2500, 5000, 7500, 10,000 ft 

V = 150, 175, 200, 225, 250, 275, 300, 325, 350 knots 

and calculate the rms percent error between the table lookup values and the analytic 
drag equation, where D is in pounds, h is the altitude in feet, and v the speed 
in knots. By means of an extremum-finding subroutine (ref. 12) the parameters are 
adjusted to minimize the rms percent error, which results in an rms error of 2.387%. 
This apparent error is primarily due to the desirable smoothing of the drag curves 
in the flap deployment region. The polynomial fit agreed within 0.8% (max) with the 
drag in the 6f = 0 region. The values of the coefficients in (A14) are: 

k x = -0. 4078xl0 4 k 3 = 0.1146 k $ = 0.3830xl0 9 

k, = 0.2429xl0 -4 k, = -0.2005xl0“ 4 k c = 0.4227xl0" 4 

Ik 6 

The resulting drag curves are shown in figure A3, which may be compared to the table 
lookup data in figure A1 . 


Fuel flow f(T,v,h) 

Reference 1 presents fuel flow data in tabular form 

i = f (Mach,T ,h) 
n n n 


where 


f 

n 

f 

n 

f 


is a normalized fuel flow rate 

= f (T/e) 
n 

= fuel flow rate 


(A15) 


(A16) 
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Figure A. 3.- Drag versus speed from algebraic equation. 
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(A17) 


T = thrust for three engines = 6 T r 
where T n is the normalized thrust for three engines. 

/e = square root of the temperature ratio 

= temperature at altitude h/temperature at zero altitude (A18) 

6 = pressure ratio = pressure at altitude h/pressure at zero altitude (A19) 


It is our task, first to determine a functional relationship for fuel flow that 
involves only T,v,h, second to determine approximate values for the parameters, and 
third to minimize the error between the analytic representation and the full table 
lookup of the fuel flow when using the standard 1962 atmosphere. 

The f n entries in the three-dimensional table are normalized in such a way that 
they are almost altitude-independent. For the first two steps of the task we are 
therefore using only the table for h = 5000 ft. In figure A4 the f n values multi- 
plied by the normalized thrust T n are plotted over the range of Mach numbers of 
interest because this quantity is the actual fuel flow normalized by atmospheric 
constants . 

From (A16) and (A17) 

T f = T f/(T/6) = f/(6/0) (A20) 

n n n 


The relationship shown in figure A4 can be approximated as 


T f 


n n 


0.45 + M + 


3.75 + 3.35 M 
30,000 n 


0.5(1 - M) 2 
30,000 / n 


(A21) 


(A20) and (A18) used in (A21) give the actual fuel flow 


r r/^Tn , .. . 3.75 + 3.25 M T 0.5(1 - M) /T\ 2 ' 

f = 6/6 [0.45 + M+ 6 + 30,0002" 


(A22) 


Fuel flow is still a function of the pressure ratio, temperature ratio, and 
Mach number, which are all functions of altitude. The altitude dependence for the 
normal atmosphere is presented in table Al. 


TABLE Al.- NORMAL ATMOSPHERE 


Altitude 

ftxlO -3 

T, 

K 

P, 

psf 

P, 

slugs /f t 3 xlO 4 

6 = P/P 0 

o 

CL 

Q. 

II 

D 

o' 1 / 2 

Sound 
speed , 
f ps 

0 

288.15 

2116 

23.77 

1.0000 

1.0000 

1.0000 

1116.9 

2 

284.19 

1968 

22.41 

.9298 

.9428 

1.0299 

1109.2 

4 

280.23 

1828 

21.11 

.8637 

.8881 

1.0611 

1101.4 

6 

276.27 

1696 

19.87 

.8014 

.8359 

1.0938 

1093.6 

8 

272.31 

1572 

18.08 

.7428 

.7860 

1.1279 

1085.7 

10 

268.35 

1455 

17.55 

.6877 

.7385 

1.1637 

1077.8 
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We can see from table A1 that the speed of sound decreases with altitude, so that the 
speed expressed in units of Mach (M) increases with altitude 

M = v(l + 0.363 x 10~ 5 h) /1116. 9 (A23) 

where v is given in feet/second true airspeed. The pressure ratio can be expressed 
as 


6 = 1/(1 + 0.45 h/10,000) - 1 - 0.31 h/10,000 


(A24) 


and the square root of the temperature ratio as 

/S = 1 - 0.363xl0" 5 h 


(A25) 


Using (A22)-(A25) allows us to obtain f when T, h, and v are given. We can 
rewrite (A22) as 

f « C 0 (h,v) + C x (h,v)T + C 2 (h,v)T 2 (A26) 

where each of the C's are functions of the atmospheric relations. When (A23) 
and (A24) are used to replace the atmospheric relations in (A22) and when all higher- 
order terms except terms involving hv are neglected, all C's end up in the same 
form: 


C. 

3 




+ c ) 
** 3+2 


+ W 1 


+ c 4 j+ 4 fc) v ; 


j = 0,1,2 


(A27 ) 


The numerical values of the twelve parameters in (A27) were used as initial 
values in the extremum-finding subroutine (ref. 12) where the parameters are adjusted 
to minimize the rms error between table lookup values for the 200 points given by all 
combinations of 


V = 150, 200, 250, 300, 350 knots 
h = 0, 2000, 4000, 6000, 8000, 10,000 ft 

while using all appropriate fuel flow tables, as well as the exact nominal atmospheric 
model. 


T = 7000, 9500, 12,000, 14,500, 17,000, 22,000, 24,500 lb 


and the values from the analytic expressions (A26) and (A27). This resulted in an 
rms error of 0.039 lb/sec. The values are: 



For idle thrust 


0.47537 
-0.24360*10"“ 
0. 17702*10" 2 
-0.28995*10"“ 


c 5 = 0. 10823*10" 3 
c 6 = -0.80509*10" 5 
c 7 = 0.15898*10" 6 
c e = -0.12439*10-“ 


c 9 = 0 . 91747*10 -9 
c 10 = 0.79644*10-“ 
c lx = -0.85521*10" 12 
c 12 = -0.45624*10"“ 


c 1 = 0.2528 
c 2 = -0.4*10"“ 
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c 3 = 0.147X10 -3 
c 4 = -0.363*10 -5 

and c 5 through c 12 are zero. These data are obtained from a table of fuel flow at 
engine idle. No data are available for an intermediate thrust of 0—7000 lb » which 
fortunately are not required for this problem. An example of the resulting fuel flow 
versus thrust function is shown in figure A5. 



0 10 20 30 


T. K-lb 

Figure A. 5.- Fuel-flow model of the 727-100. 
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APPENDIX B 


In further development we need certain time and other derivations which we 
enumerate here. Let's assume that T = constant and u = time variable. 

F = [c^l + c 2 h) + c 3 (l + c 4 h)v] + [c 5 (l + c 6 h) + c 7 (l + c e h)v]T 
+ [c 9 (l + c 10 h) + Cll (l + c 12 h)v]T 2 
F v = c 3 (l + c 4 h) + c 7 (l + c 8 h)T + c 13 (l + c 12 h)T 2 

K = ^ C 3 C *» + C 7 c 0^ + C 11 C 12^- ^ = F vh^ 

F h = (c 1 c 2 + c 3 c 4 v) + (c 5 c 6 + c 7 c 8 v)T + (c 9 c 10 + c^c^v)! 2 
^h = 3 C 4 + C 7 C 8 T + C 11 C 12 T ^ = F vh^ 

D = k 3 (l + k 2 h) + k 3 (1 + k 4 h)v 2 + k 5 (1 + k 6 h)/v 2 

+ [k 3 (1 + k 2 h) /2 + k 5 (1 + k 6 h)/v 2 ]u 2 = D x + D 2 u 2 

D = 2D,u 
u 2 

D v = 2k 3 (l + k 4 h)v - 2k s (l + k 6 h)(l + u 2 )/v 3 = D lv + D 2 v u 2 

D v = 2k 3 (l + k 4 h)v + 2k 3 k 4 hv + 6k 5 (l + k 6 h)(l + u 2 )v/v 4 

- 4k 5 (l + k 6 h)uu/v 3 - 2k 5 k 6 h(l + u 2 )/v 3 

= k 3 k 2 + k 3 k 4 v 2 + k 5 k 6 /v 2 + (k 3 k 2 /2 + k 5 k 6 /v 2 )u 2 = D lh + D^u 2 

= 2k 3 k 4 vv + (k 3 k 2 + 2k 5 k 6 /v 2 )uu - 2k 5 k 6 (l + u 2 )/v 3 v 

D = D h + D u + D v 
h u v 



u = 0 ; u = u m 


shall 

(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 
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APPENDIX C 


• • 

For the minimum principle to apply, the velocity set (h,i^,v,F) must be convex 
(refs. 13 and 14). Using equations (3) and (5) for our particular example, fig- 
ure Cl (a) shows that for descent, the velocity set is not convex in the region of 
zero thrust, due to increased drag at nonzero bank angles. In this case a relaxed 
(chattering) controller may result in better performance if maximum deceleration will 
improve performance. To convexize the velocity set we allow the possibility of a 
chattering bank angle and define a combined control. 

u = au - (1 - a)u = u (2a - 1) (Cl) 

s s s 

where u s is the magnitude of the tangent of the bank angle u s = | tan <f> | and a 
is the chatter parameter (fig. Cl(b)). 

a = per unit time the bank angle spends at u g 
1 - a = per unit time the bank angle spends at -u 
then (4) changes to 


= -gu g (2a - l)/v 


(C2) 


and (12) becomes (assuming n = 0, off the speed limit) 

H = F - Cv cos <p + X^vy - X^gu g (2a - l)/v + X^g(T - D x - D 2 u^ - Wy)/W (C3) 


Note that the drag is only a function of u| and does not depend on the value of 
the chatter parameter a. Now 


i* = 


if 

if 


X, < 0 

X > 0 
'll 


(C4) 


and (20) and (21) become 


u* = 
s 


m 


if 

A 

and 

X 

V 

< 0 

(C5) 

if 

V £ u„ 

m 

and 

X 

V 

< 0 



v = [ (2a - l)WX^]/[2X v D 2 v] 


(C6) 


We note that when a is nonsingular (49)-(51) give the same result as (21). In 
this case, (C2) and (C3) revert to the original equations (4) and (12) and we can 
solve the optimal descent problem as if the original velocity set were convex. 
However, a is linear and we must consider the singular-a solution. 


H = -2gX.u /v = 0 ->■ X = 0 
a ° ^ s ^ 


( 07 ) 


39 





(C8) 


(H ) = 2gCu sin = 0 -*■ \p = 0 
Ot s 

GO = -2g 2 u 2 C cos ij»(2a - 1) = 0 a = 1/2 (C9) 

cx s 

/ u \ / 2 2 P , / ^ n Since C > 0 \p = 0 

(-H ) = 4u g C cos \p/v > 0 r (CIO) 

m cx o , 

but not \p = ir 

But from (20) and (C5), u s = 0 when A^ = 0. This is equivalent to the nonchatter- 
ing straight descent and again is included in the earlier solutions, which do not 
require a chattering bank angle. 

A different approach to make the velocity set convex is to apply appropriate 
constraints. This will cut off portions of the earlier velocity set. A reduced 
velocity set in general means reduced performance. We will therefore reduce the 
velocity set in such a way as to minimize the excised area. To do this we impose an 
acceleration constraint, which limits deceleration to a constant value for each (v,h) 
such that it equals the maximum value which can be obtained at zero bank angle 
(fig. Cl(c)). The resulting deceleration constraint is a mixed state and control 
constraint 


v > -g[D x (v,h) + W Ymax ]/W (Cl 1) 

In the following development we only express the changes from the original set of 
equations, where we express the original equations as before and star (*) the new 
values. Then from (12) and (C 1 1 ) 

H* = H + y{-g[D (v,h) + Wy ] /W - v} (C12) 

-J- max 


Inserting the equation for v and simplifying gives 

H* = H + u(-g[T - D 2 (v,h)u 2 - W( Y - Y n/W) (C13) 

III 3.X 

the adjoint equations (13), (14), and (15) become 


= A h - 

ygD^u 2 /W 

(C14) 

X* = X, 


(C15) 

X* = X - 

V V 

ygD V2 U 2 /W 

(C16) 


On the deceleration bound we have T - D ? (v,h)u 2 - W(y - y ) = 0 such that 

2 1 ’max 

Y = Y - D 9 u 2 /W + T/W (Cl 7) 

max 2 


We are interested in the decelerating descent at T = 0, then y is just decreased 
sufficiently to compensate for the additional drag D 2 u due to a nonzero bank angle. 
From H* = 0 = Hy + yg we have 
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1 A„ - X, v/ g >0 on the deceleration constraint 

v h (CIS) 

0 off the deceleration constraint 

The >0 requirement for y on the acceleration constraint is a necessary condition 
for optimality. Thrust control 


H* = H x - yg/W 


Hence the only change from (17) is 

t = -[C 1 (h,v) + (X v - y)g/W] /[2C 2 (h,v)] 


(C19) 


and for bank angle control 


H* = H u + y2gD 2 (h,v)u/W 


(21) becomes 


v = -A^W/[2(A v - y)D 2 (h,v)v] (C20) 

Note in (C19) and in (C20) A v - y can be replaced via (C18) by 

A v - y = A h v/g (C21) 

We are now in a position to compare the solutions for the convex set with those of 
the "natural 11 set. If the latter solutions have a better performance we shall accept 
them, since the acceleration constraint was only artificially introduced to satisfy 
the minimum principle rather than representing passenger or equipment requirements. 

The above constraints, which only apply to the descent, did not affect the performance 
(less than 0.2 lb worse than the unconstrained optimal). 
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